Graphite intercalation compound KCs revisited: a Icey to 
graphene 



A. Gruneis\ C. Attaccalite^, A. Rubio^, D. Vyalikh^, S.L. Molodtsov^ J. Fink^ R. FoUath^, 
W. Eberhardt^, B. Buchner\ T. Pichler^ 

^IFW Dresden, P.O. Box 270116, D-01171 Dresden, Germany 

"^Nano-Bio Spectroscopy Group and European Theoretical Spectroscopy Facility (ETSF), Departa- 
mento de Fisica de Materiales, Unidad de Materiales Centra Mixto CSIC-UPV/EHU, Universidad 
del Pais Vasco, Avd. Tolosa 72, E-2001 8 Donostia, Spain 

^Institut far Festkorperphysik, TU Dresden, Mommsenstrasse 13, D-0 1069 Dresden, Germany 

^BESSYn, Albert-Einstein-Str 15, D- 12489 Berlin, Germany 

^Faculty of Physics, University of Vienna, Boltzmanngasse 5, A- 1090 Vienna, Austria 

Electrons in isolated graphene layers are a two-dimensional gas of massless Dirac Fermions. 
In realistic devices, however, the electronic properties are modified by elastic deformations, 
interlayer coupling and substrate interaction. Here we unravel the electronic structure of 
doped graphene, revisiting the stage one graphite intercalation compound KCs using angle- 
resolved photoemission spectroscopy and ab-initio calculations. The full experimental dis- 
persion is in excellent agreement to calculations of doped graphene once electron correlations 
are included on the GW level. This highlights that KCs has negligible interlayer coupling. 
Therefore Dirac Fermion behaviour is preserved and we directly determine the full experi- 
mental Dirac cone of doped graphene. In addition we prove that superconductivity in KCs is 
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mediated by electron-phonon coupling to an iTO phonon, yielding a strong kink in the quasi- 
particle dispersion at 166 meV. These results are key for understanding, both, the unique 
electronic properties of graphene and superconductivity in KCg. 

The recent discovery of two-dimensional meta-stable graphene sheets has sparked enormous 
interest in their low-energy electronic structure Graphene samples are in general obtained by 
three methods: (1) repeated peeling of a graphite single crystal ^, (2) growth by chemical vapour 
deposition on Ni(lll) ^'^ and (3) precipitating few-layer graphene from SiC ^. Angle-resolved 
photoemission spectroscopy (ARPES) has been proven to be a key tool to determine the electronic 
structure of one and few layer graphene ^ and graphite "^'^'1°. A major problem for the investigation 
of substrate based graphene is that there is a significant modification of the electronic structure due 
to interaction with the substrates yielding charge transfer and hybridisation ^. 

One way to overcome this problem is to measure single crystalline graphite, which has no 
substrate interaction. In this case however, a kz dispersion of two tt valence bands ^ and a small 
gap were observed because of the AB stacking sequence of graphene layers. Another important 
issue in both systems is the renormaUzation of the bare electronic band structure due to doping de- 
pendent electron-electron correlation ^, electron-phonon coupling (EPC) or electron-plasmon 
coupling 

In order to circumvent the problems of substrate interaction, strong bilayer splitting and 
electron-electron correlation we revisited stage one graphite intercalation compound (GIC) KCg. 
GICs have been at the focus of intense research in the last four decades because they have a wide 
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range of tunable electronic properties . Especially, for stage I GICs superconductivity was ob- 
served with transition temperatures Tc ranging from below 1 K for alkali metal intercalation (e.g. 
0.55 K for KCg up to 11.5 K for CaCe In both cases EPC is the superconducting pairing 
mechanism In stage I alkali GICs the graphene layers have AA stacking and only one tt 
conduction(valence) band and can thus be considered as a doped graphene layer sandwiched in be- 
tween two positively charged plates (Supplementary Fig. SI). Thus the low-energy band structure 
of, both, stage I GICs and graphene are described by a 2x2 tight-binding (TB) Hamiltonian ^^'^^ 
resulting in a linear tt band dispersion close to the crossing point of the valence and conduction 
bands. The electronic band structure of GICs was also calculated by density functional theory in 
the local density approximation (LDA) Experimentally, preliminary studies on ARPES of GICs 
were reported Until now the details in the low energy quasiparticle (QP) dispersion of GICs 
regarding the superconducting coupling mechanism are not identified. Furthermore the issue of 
whether the charge transfer to graphite is complete or partial 22,23,28 never resolved. 

In this work we revisited the electronic structure of KCs GICs using a combination of ARPES 
and ab-initio calculations. We proof a complete charge transfer and find excellent agreement to 
ab-initio GW calculations including electron-electron correlation. This highlight that a rigid band 
shift model of graphene is applicable for KCg and electron-electron correlation play a crucial role 
in the QP dispersion. Hence we unravel for the first time the full experimental Dirac cone of 
graphene without modifications by substrate interactions and directly determine the corresponding 
momentum dependent Fermi velocites {vp) in the valence and conduction bands. Thus our results 
are key for understanding the unique doping dependent electronic transport properties of gated 

3 



graphene layers in nanoelectronic devices. In addition we show a detailed study of the direction 
dependent renormalization of the QP dispersion at low binding energy that highlights the coupling 
to an iTO phonon at 166 meV. We directly show that the coupling to this phonon, which is also 
responsible for the double-resonance Raman process, is the major contribution to the supercon- 
ducting pairing in KCs- 

The KCs crystal structure is given by individual graphene sheets separated by layers of potas- 
sium as shown in Fig.[T]^a). This compound was synthesized in-situ by evaporation of potassium as 
described in detail in the Methods section. The fully intercalated graphite crystals have a character- 
istic golden colour as shown in Fig.[TJb). We now carry out a detailed analysis of the QP dispersion 
of KCg as measured by ARPES. In Fig.[T]^c) we show energy dispersion curves (EDC) and cuts of 
the TT bands that intersect the comers of the Brillouin zone (BZ). The cuts are done close to the ky 
direction [see coordinate system in Fig. Qb)]. We also measured the photon energy dependence 
and, beside changes in the matrix elements, do not observe a dispersion perpendicular to the lay- 
ers. Fig.[T]^d) shows equi-energy cuts of the raw photoemission intensity. When moving from the 
cut dX Ep io lower energies the circumference of the equi-energy contour of the conduction band 
becomes smaller. Evidently we have one vr valence and one tt conduction band and they meet in 
one point. Therefore in akin analogy to graphene the crossing point at 1.35 eV is from now on 
called the "Dirac point". Decreasing the energy of the cut below 1.35 eV enlarges the equi-energy 
contour of the valence band. Comparing two contours with the same distance from the Dirac point 
[e.g. cuts at Ep and 2.7 eV in Fig. [TJd)], it can be seen that the tt valence band is steeper and 
has less trigonal warping. Such a result is related to the fact that the self-energy corrections to the 
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trigonal warping effect increase with doping ^. Interestingly, as depicted in Fig.yjd), the photoe- 
mission matrix elements of the valence and conduction bands are complementary yielding a high 
photoemission intensity in 1/3 and 2/3 of the BZ, respectively. From the contour at Ep in Fig. ^d) 
we determine the number of carriers rie by integrating the volume inside the Fermi surface. This 
yields Ue = 7.2 x 10^^ electrons cm^'^, close to a full charge transfer (Ue = 8.9 x 10^^ electrons 
cm^^). The remaining deviations can be explained by a slightly lower stochiometry, i.e. Kg.ssCg. 
From these results, concomitant with the absence of any Fermi surface of K 4s states close to the 
r point (Supplementary Fig. S7) we can safely state a complete charge transfer from potassium in 
agreement with previous results 24-^7 ggg^ Fig.|2]the ARPES intensity maxima at Ep 

and the linear dispersion close to the Dirac point are very well described by a TB fit including third 
nearest neighbours (NN) (Supplementary Table 1). First NN TB calculations are inadequate and 
can not describe the experimental equi-energy contours. This new TB fit reproduces very well the 
experiments shown in Fig. |2] and will be also used for the self-energy analysis of the EPC to be 
described below. 

In order to further understand the underlying band structure on an ab-initio level we com- 
pared the experimental results to calculations at six different levels comprising LDA and GW ^""^^ 
calculations of (un)doped graphene and KCg. As seen in the Supplementary in Fig. S2-S6, band 
structures calculated at the level of LDA are too flat yielding an incomplete charge transfer with 
a partially occupied K 4s band in the vicinity of the T point. Calculations on the GW level in- 
cluding full and partial self consistency in G are in much better agreement to the experiment. The 
obtained results for the QP dispersion at K are identical but only the latter explains the lack of pho- 
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toemission intensity at the T point in KCg (Supplementary Figs. S7-S8). As depicted in Fig. ^a) 
the observed linear dispersion is reminiscent of the underlying structure of the graphene parent 
compound and in perfect agreement to the GW calculations of graphene. This also holds for the 
equi-energy contour of the photoemission intensity at Ep shown in Fig.|2];b). These results clearly 
indicate that electron-electron correlation is crucial to explain the size of the Fermi surface and the 
observed trigonal warping. I.e. electron-electron correlation is at the heart of the band structure 
the underlying (doped) graphene layers. Hence our results unambiguously highlight the following 
facts inherent to the electronic structure of doped graphene without substrate interaction: (i) the 
shape of the Dirac Fermions depends on the doping level and the equi-energy contours above and 
below the Dirac point are anisotropic; (ii) the trigonal warping increases as doping increases; (iii) 
in KCg each potassium atom transfers one electron to the graphene layers. 

Most importantly let us now turn to a quantitative assessment of the QP bandstructure around 
the Dirac point, i.e. unravelling for the first time the full experimental Dirac cone. For this pur- 
pose, we extracted the photoemission maxima for binding energies between Ep and 3 eV in steps 
of 10 meV. The results are depicted in Fig.[3];a). It is obvious that the two-dimensional band struc- 
ture close to the Dirac point is linear and two bands meet each other in one point. This highlights 
that KCs indeed consists of doped sheets of graphene and can be well described with the aforemen- 
tioned TB and GW approaches. We point out this this is in contrast to all previous experiments 
on substrate based graphene, clearly observing the opening of a gap in the electronic and structure 
of epitaxial graphene on Ni(l 1 1) ^ and SiC Although the origin of the gap in graphene on SiC 
is under debate it existence is accepted. The opening of a gap causes a breakdown of Dirac 
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Fermions which has profound limitation on the observation of relativistic physics in graphene. 
Therefore the spectroscopic investigation of doped graphene layers in KCg provides an elegant 
solution to this problem. 

Hence, by analyzing the full two-dimensional dispersion of KCg around the Dirac point we 
can thus actually learn about the physics of graphene. In order to further underline the validity of 
this analogy we analyze the experimental direction dependent in the region of 0.1 eV above 
and below the Dirac point. These results for the valence and conduction band along with the GW 
calculations of doped graphene are depicted in Fig. ^h) and Fig. ^c), respectively. The vp has a 
maximum(mininimum) in KT(KM) direction. The very good quantitative agreement highlights 
the validity of our approach. In close similarity to our results on pristine graphite ^, we observed 
that electron-electron correlation is crucial to explain the band structure of doped graphene layers 
in KCg. 

We now turn to the analysis of QP bandstructure of KCg close toEp and their implications 
for superconductivity. In addition to the linear dispersion, as can be see in Fig. [TJ^c) Fig. ^a) and 
Fig. [3]; a), the vr conduction band of KCg is strongly kinked at 166 meV. This is in agreement to 
previous results on doped graphene By an accurate comparision to the graphite phonon disper- 
sion relation ^^"^^ we can unambiguously assign the kink to a coupling to the in-plane transverse 
optical (iTO) phonon branch near the K point (see the Supplementary Table 2 for the calculated 
phonon frequencies including non-adiabatic effects). This agrees perfectly with both the energy of 
the kink and also with the facts that the EPC matrix element with the iTO phonons is much stronger 
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than with T point phonons and that the phonon density of states is strongly peaked due to the 
flat dispersion of the iTO branch around K. The photohole decay process is shown schematically 
in Fig. |4];a). It is clear that a photohole can relax to a lower binding energy state in two ways by 
intravalley and intervalley scattering. For the case of intravalley scattering the photohole scatters 
around K or K' points and relaxes by emission of T point phonons. However, as their energy 
does not fit with the measured kink we disregard this mechanism. Therefore the relevant process is 
related to intervalley scattering where the photohole scatters between K and K' points and relaxes 
by emitting phonons close to the K point. In Fig. ^h) we illustrate this mechanism in the 2D 
BZ of graphene with the exchange of a phonon with wavevector qp^. Interestingly, a very similar 
process for scattering of photoexcited electrons is responsible for the D and G' band in the double 
resonance Raman process in sp^ hybridized carbon materials. 

Once we have identified the phonon mode responsible for the observed kink it is important to 
accurately determine A, the EPC constant In recent works A was estimated from the change 
in the slope of linear bare electronic bands; however this is not justified as the proper bare band 
structure at the Ep has to be considered, as was clearly pointed out in a recent theoretical work 
In the present work we avoid this problem by considering the total self-energy S(-E) from which A 
can be directly obtained by a simple energy derivative at Ep A = — ^-^'^^^(-^)) | ^^^^ . To do so we fit 
the momentum distribution curves (MDCs) of the QP band structure by Lorentzians. The positions 
of the maxima and the widths of the MDCs are directly proportional to the real and imaginary part 
of S, i.e. Re(S) and Im(S), respectively (see Supplementary Figs. 9-10). The QP dispersion in the 
energy region close to the kink is shown in detail in Fig. Qc) along with the measured maxima of 
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MDCs and the bare band dispersion (the cut is in ky direction in between KM and KT and thus 
corresponds to an averaged A). In the simple case of coupling to one Einstein phonon Re(S) has a 
peak and Im(S) has a jump at the energy of the coupling phonon mode Indeed this is what we 
observe in the Figs. |4];d,e) due to the strong coupling to the iTO phonon from K point at 166meV. 
Here we emphasize that the evaluation of A from Im(S) does not depend on the details of the 
dispersion of the bare bands. For the analysis of A from Re(S) we used the TB calculation of the 
bare bands which gives the correct curvature in KM and KT directions. Thus our analysis does not 
suffer from an overestimation of A in KM direction. Moreover, from the slope of Re(S) and from 
the height of the jump in Im(E) it is possible to evaluate the EPC A^^: the values we got are A = 0.4 
and A = 0.3 from the analysis of Re(S) and Im(S), respectively. This procedure can be used to 
obtain the wavevector dependence of A^. In Fig.^f) we show that the EPC obtained from both the 
Re(S) (squares) and Im(S) (circles): we thus obtain that A^ is maximum close to KM direction 
and minimum in the KT direction. As can be seen in the Fig.^f) the resulting A from Re(Il) and 
Im(E) are in good agreement regarding the anisotropy of the coupling although the actual size of 
the coupling differs by about 20 %. The minor anisotropy obtained in previous DFT calculations 
compared to the one extracted from the present work [Fig,|4];f)] and other experimental works 
raises some concerns related to the need of going beyond LDA to describe exchange and correlation 
effects in the QP energies and wavef unctions. As we stated above, self-energy corrections at the 
GW level are needed to describe the observed absence of ARPES intensity at T, in contrast to 
what is predicted by LDA. Another possible explanation to partly explain the strong anisotropy in 
Afc observed experimentally might be the phonon trigonal warping effect Thus further work is 



9 



needed to address the phonon dispersion relation and the impact of the GW self-energy corrections 
on EPC which we are currently working on. We have also analyzed /3, i.e. the coefficient 
of the quadratic energy dependence of Im(Il) which is depicted as crosses in Fig. ^f) (see also 
Supplement Fig. 10). Interestingly, A and /? have the same anisotropy, which reflects the anisotropy 
of the QP bandstructure at Ep. 

The strong EPC limits the maximum bias current through graphene based nanoelectronic 
devices and is responsible for the strong double-resonance Raman scattering in sp^ bonded carbon 
materials. In classical BCS theory, with a pairing mechanism based on EPC it is also one of the 
limiting factors for the transition temperature of KCg. The evaluated transition temperature using 
the McMillan formula '^^ with the high phonon temperature 6 = 1926 K and the average EPC 
constant of A = 0.45 [Fig.^f)] and a screened pseudopotential /i* ~ 0.14 would allow a Tc up 
to 6 K. These remaining differences to the experimental Tc can be explained by either uncertainties 
in /i* or by additional not yet known pair breaking effects. 

In conclusion we have synthesized KCs in-situ and performed ARPES measurements. We 
have found a complete charge transfer and that the QP dispersion of KCg is strongly modified 
in the presence of EPC to the iTO phonon branch at K. The value of the coupling constant is in 
good agreement to the experimental value of Tc and hence determines the superconducting pairing 
mechanism. Interestingly, this iTO phonon branch is also responsible for the strong double reso- 
nance Raman scattering in sp^ bonded graphite materials. The QP band structure of KCg is in good 
agreement to ab-initio calculations for graphene after including self-energy corrections on a GW 
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level. We thus conclude that the conical band dispersion of KCg closely resembles that of Dirac 
Fermions of graphene after applying a gate voltage. Our results circumvent the problems which 
are associated with strong substrate interaction of graphene ^'^^'^^ or with interlayer interaction ^, 
both of which cause a breakdown of Dirac Fermions. Most importantly, we unravelled for the first 
time the full experimental Dirac cone of doped graphene and directly determine the corresponding 
momentum dependent vp va^ the valence and conduction band as such providing crucial input for 
the understanding of the unique electronic and transport properties of graphene. 

Methods 

Experiments were done at BESSY n using the UE112-PGM2 beamline and a Scienta RS 4000 
analyzer yielding a total energy resolution of 15 meV and a momentum resolution better than 
0.01 A~^. Natural graphite single crystal samples were mounted on a three axis manipulator and 
cleaved in-situ. Intercalation was performed with the sample at room temperature by evaporating 
potassium metal from a commercially available SAES getter source. In order to check the dop- 
ing level, we measured ARPES after each intercalation step. The proof that we reached stage I 
was given by the appearance of only one tt valence band [instead of 2(3) valence bands for stages 
n(in)]. In addition stage I compounds are identified by their characteristic golden color. After full 
intercalation, KCg was immediately cooled down by liquid He to 25 K and measured. The calcu- 
lations of the electronic dispersion of graphene in a slab geometry (d = 20 a.u.) are performed on 
two levels. First, we calculate the Kohn-Sham band-structure within the LDA to density-functional 
theory Wave functions are expanded in plane waves with an energy cutoff at 25 Ha. Core elec- 
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Irons are accounted for by Trouiller-Martins pseudopotentials. In the second step, we use the 
GW approximation ^^'^^ to calculate the self -energy corrections to the LDA dispersion For the 
calculation of the dielectric function e(a;, q) we use a Monkhorst-Pack k grid sampling 36x36x1 
points and bands up to 70eV (namely 50 bands) of the first BZ. In a further step we performed GW 
calculations including partial self consistency in G. Note, that this is a conserving approximation 
whereas the single shot GqWo is not. 
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Figure 1 (a) Crystal structure of KCs- (b) Photo of intercalated KCs single crystals with 
golden colour, (c) As measured EDC cuts through the corner of the 3D BZ. (d) Equi- 
energy contours of the photoemission intensity taken at 48 eV photon energy for four 
binding energies. 

Figure 2 (a) ARPES scan measured close to the ky direction along with the TB fit of the 
bare-band dispersion (black) and the GW calculation for electrostatically doped graphene 
(blue) and undoped graphene (green circles), (b) Symmetrized equi-energy contour for 
E=0.24 eV and maxima (crosses) along with the TB fit and the GW ab-initio calculations. 

Figure 3 (a) Experimental Dirac cone from the observed photoemission intensity max- 
ima (denoted as dots). Measured and calculated {GW) values of vp for (b) electrons and 
(c) holes around the Dirac point. 

Figure 4 (a) Schematics of the photohole scattering process from K to K' by EPC to 
a K point phonon with energy huph{qph)- The initial and final photohole states around K 
and K' are denoted by circles. The arrow denotes the exchanged phonon with wavevector 
qph and energy %ojph{qph) which is emitted by photohole relaxation, (b) 2D BZ indicating 
the phonon wavevector qph. (c) measured ARPES intensity in the region of the kink. A 
straight dashed line denotes the bare band dispersion and a kinked dotted line are the 
experimental MDC maxima, (d) real and (e) imaginary part of the self-energy, (f) the 
angular dependence of the EPC constant A and the electron-electron correlation constant 
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(5; black squares(red circles) denote A extracted from the measured real(imaginary) part 
of the self-energy. j3 is denoted by blue crosses. 



24 



